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Abstract 



We discuss the discovery reach of LEP2 for the Higgs sector of a general extension 
of the MSSM including a single gauge singlet field. This change introduces a new 
quartic Higgs boson self-coupling which can increase the masses of the CP-even states, 
and also allows mixing between singlet and non-singlet states which can reduce the 
couplings of the mass eigenstates to the Z. The lightest CP-even Higgs boson is 
bounded by a parameter A which takes a maximum value A max « 136 — 146 GeV for 
top mass 150 — 195 GeV. We generalise the discussion of the bound to include the 
entire CP-even spectrum and show how experiment may exclude values of A smaller 
than some A min . CP-even Higgs boson searches at LEP2 will be able to exclude 



A m j„ w 81 — 105 GeV, depending on the machine parameters. We also present 
exclusion plots in the uia — tan (3 plane, based on an analysis of CP-even, CP-odd 
and charged Higgs production processes at LEP2. 
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1 Introduction 



The minimal supersymmetric standard model (MSSM) Q has a Higgs sector consist- 
ing of two doublets if i and H2 coupling to the down-type quarks and charged leptons, 
and to the up-type quarks respectively. The particle content, other than Goldstone 
bosons, is then two CP-even states h, H; one CP-odd state A; and a charged scalar 
H^. The requirements of supersymmetry and gauge invariance then constrain the 
quartic Higgs coupling in terms of the gauge couplings, and we are left with only 
two free parameters to describe the whole Higgs sector. These are conventionally 
taken to be m^, the mass of the CP-odd state, and tan/? = Vijv\ where ^ = (H®), 
and v\ + v\ = v 2 = (174GeV) 2 . It is then straightforward to derive the formulae 
m\ < Mz and m#± > My/ where the first of these in particular is greatly affected by 
radiative corrections H. Thus the MSSM is quite constrained and leads to the usual 
LEP2 Higgs discovery limits in the iua — tan /3 plane. 

If the Higgs sector of the MSSM is extended by the addition of further particle 
content then much of this predictivity is lost. This occurs for two reasons. Firstly, the 
constraint on the quartic couplings in the Higgs sector is not purely a result of super- 
symmetry but also an artefact of having only doublets in the Higgs sector, making it 
impossible to introduce extra Yukawa couplings through the Higgs superpotential in 
a way consistent with gauge invariance. Increasing the particle content and including 
such Yukawa couplings destroys both the upper bound on and the lower bound 
on m#±. Secondly, the extra states which we can introduce will mix with the states 
present in the MSSM, and can alter both their couplings and their masses. 

For example, in the next-to-MSSM (NMSSM) where there is an extra gauge singlet 
state N which only couples through the Higgs superpotential, and where there are 
only trilinear terms in the superpotential 0, there are now three CP-even neutral 
Higgs scalars and two CP-odd neutral Higgs scalars, due to the real and imaginary 
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components of the additional singlet scalar. In this model the lightest CP-even neutral 
scalar may be significantly heavier than in the MSSM, and to make matters worse 
this lightest CP-even scalar may have diluted couplings to the Z boson due to the 
admixture of singlet component ffl . However, it is possible to derive a bound on the 
lightest CP-even state, both in this and in more general models ||, ^, f|. Furthermore, 
in the limit that the lightest CP-even Higgs boson is completely decoupled (and hence 
we might think that the bound does not tell us anything about states which could 
be detected), the bound applies instead to the second lightest CP-even Higgs boson, 
while for light states which are merely weakly coupled one can derive precise bounds 
on their heavier partners. These become closer to the bound on the lightest as the 
singlet component of the lightest becomes greater [/J. 

Our intention in this paper is to discuss how much of the parameter space of such 
a general SUSY model can be covered by searches at LEP2, comparing them closely 
with the corresponding results for the MSSM. We shall consider the most general 
possible model with a gauge singlet and show that LEP2 can discover CP-even Higgs 
states in any part of the tra — tan (3 plane, and can exclude significant parts of this 
plane even with the most general possible mixing. 

We begin with a short review of the MSSM and its simplest extensions and how 
they affect the mass matrices. In section 4, we discuss the conventional bound on 
the lightest CP-even state in this model. In section 5, a number of search strategies 
are discussed, and we consider how effective the various MSSM searches are in an 
extended model. Finally, in section 6 we present results applicable to LEP2. Section 
7 is the conclusion. 
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2 The MSSM 



Here the superpotential is: 

Wmssm = -^H 1 H 2 + ... (2.1) 

where #i# 2 = H^H® - H{H^. 

The superpotential leads to the tree-level Higgs potential: 

Vmssm = rnllH.l 2 + m 2 2 \H 2 \ 2 + m 2 12 (H 1 H 2 + H.c.) 

+ l( g 2 + g ' 2 )(\H 1 \ 2 -\H 2 \ 2 ) + ^\H* 1 H 2 \ 2 (2.2) 

where m 2 = m 2 Hl + /i 2 , m 2 = m 2 H2 + fi 2 , and m\ 2 = —Bjj, where B and m 2 H . are the 
usual soft SUSY-breaking terms. 

The tree-level CP-odd mass-squared matrix is 

M\=( m H? ™\ ) (2.3) 
^ m\ 2 m 2 2 /tp ) 

M\ is diagonalised by UM A W = diag(m 2 1 , 0) where U is given by 

U=( Sp (2.4) 

and rriA is given by m 2 A = m\ 2 j (spcp) where sp = sin (3, cp = cos (5 and tp = tan (3. 

The tree-level charged Higgs mass squared matrix is 

M 2 _ ( M^s 2 + m\ 2 tp M^spcp + m\ 2 \ 
H± ~ { M 2 w s p cp + m\ 2 M 2 w c 2 + m 2 12 /tp ) 

Mjj± is diagonalised by UM^±W = diag(m^±,0) where U is as before and m H ± is 
given by m 2 H± = + m 2 A . 

The tree-level CP-even Higgs mass-squared matrix is, after substituting the tree- 
level expression for m\ 2 in terms of m\, 

I M|c| + m\sl -(Ml+m\)s f ,\ 



M 2 is diagonalised by VM 2 V^ = diag(m 2 H , ml) where 

y = f C a S a \ ^ ^ 

where — n/2 < a < and m 2 h1 m\ are given by 

m\ h = \[m\ + Ml) ± l -^m 2 A + M 2 ) 2 -Am 2 A M 2 c 2 w (2.8) 

and the two angles are related by: 

tan 2a = tan 2 ^™} + M f? (2.9) 

In the above conventions, the ZZh coupling has an additional factor of sin(/3 — a) 
relative to the standard model, and the ZZH coupling has a factor of cos(/9 — a). The 
ZhA coupling is proportional to cos(/9 — a). 

It is interesting to observe that if we do not diagonalise M 2 but instead rotate 
to a primed basis where the second Higgs doublet does not have any VEV, then the 
11 element of the matrix no longer contains any dependence. Since the smallest 
eigenvalue of a real symmetric matrix is less than or equal to the smallest diagonal 
term, this gives a useful upper bound on the lightest CP-even Higgs mass. This 
rotation is given by WM 2 W^ = M' 2 where 

W — ( 0/3 S ?) (2.10) 

and the bound is then given by 

m 2 h < M[l = Mf cos 2 2(3. (2.11) 
Starting from the primed basis the CP-even Higgs matrix is diagonalised by 

XM' 2 X^ = diag(m^, m 2 ) (2.12) 

where X = VWl 

Note that in this basis the additional factor in the ZZh coupling relative to the 
standard model of sm((3 — a) is simply X 2 \, and similarly the ZZH factor of cos(/3 — a) 
is just Xu. 



3 The MSSM plus a Singlet 



Now let us consider a very general Higgs superpotential of the form 

W = Wmssm + XNHiH 2 + /(AT) (3.1) 

where / is an arbitrary holomorphic function of N, the singlet field, and Hi and H2 
are the usual Higgs doublets coupling to down and up quarks respectively. Since the 
singlet does not have gauge couplings we may then write 

V = Vmssm + X 2 \HiH 2 \ 2 

+ (\(N + N)fi + A 2 |Af ) (\Hi\ 2 + \H 2 \ 2 ) 

+ (Aj^ - XNA X )H 1 H 2 + h.c. + ■■■ (3.2) 

where the ellipsis indicates terms with no dependence on Hi or H 2 . Given that Vmssm 
includes terms /j, 2 (\Hi\ 2 + l-f^l 2 ) and BfiHiH 2 it is then clear that we can account 
for the effects of all the singlet dependent terms on the upper 2x2 block of the mass 
matrices (CP-odd, CP-even, and charged) purely by redefinitions of B — ► B' and 
\i — > // in terms of the VEV of the singlet < N >= x and by including the effects of 
the A dependent quartic term. 

We further note that we could obtain identical results with an arbitrary number 
of singlets, since we can always rotate so that only one of them couples to HiH 2 and 
the remainder simply complicate the form of / and of the soft terms, which we are 
essentially regarding as arbitrary. With one singlet, the most general superpotential 
is 

k 1 

W = -nHiHz + XNH1H2 - -N 3 + -/i'N 2 + fi"N (3.3) 

3 2 

In the limit that fi,fi'fi" — > 0, the above superpotential reduces to that of the 
NMSSM0, while if N is removed it reduces to that of the MSSM. 

The mass matrix for the CP-odd scalars in the basis (Hi, H 2 , N) is a complicated 
3-dimensional generalisation of M\ in the MSSM. With only minimal risk of ambiguity 
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we shall use the same notation for such generalisations as for the MSSM. The tree-level 
CP-odd mass squared matrix becomes : 



12 
12 



V • • • 



(3.4) 



where mf 2 = —B'fj,'. The CP-odd matrix has its Goldstone modes isolated by 




(3.5) 



where U is given by 



/ Cp Sp \ 

U= sp op (3.6) 

V o o i J 

Clearly m' A is the analogue of m A of the MSSM, and w!\ = mf 2 /(spcp). The entries 
in these matrices represented by dots are complicated functions of soft terms and 
parameters from / and the extended soft potential, and since / is arbitrary they are 
essentially unconstrained. We define the matrix Va which diagonalises M\ as follows 



(V A U)M 2 A (V A Uy = diag(0,r< ,m 2 M ) 



(3.7) 



where we have taken m Al < m A2 and Va can be represented in terms of one mixing 
angle 7 



V A 



1 








\ 


















J 



(3.8) 



As regards the charged Higgs sector, the singlets obviously cannot mix with 
charged scalars, and we find that (at tree-level) the mass of the charged Higgs is 
given by 

m 2 H± = m'l + Mw — AV (3.9) 

We can immediately obtain the MSSM results by simply setting A = and removing 
all of the singlet terms from the mass matrices. Clearly a non-zero A tends to reduce 
the charged scalar masses which can be arbitrarily small. 



The CP-even mass squared matrix is now a 3x3 matrix in the basis (Hi, if 2, N) 
which may be written as 

/ M 2 z c 2 +m' 2 s 2 -{Ml + m'l - 2AV) W . \ 

M 2 =\ -(M 2 z + m /2 -2\ 2 v 2 )spCp M 2 z s} + m%<% . (3.10) 

where, as before, the dots indicate complicated entries. As in the MSSM the CP-even 
Higgs matrix is diagonalised by 

VMW = diag(m 2 1? m 2 h2 ,m 2 h3 ) (3.11) 

where V is some complicated 3x3 matrix. Note that we have now chosen to order 
the mass eigenstates as 

m hl < m h2 < m h3 (3.12) 

so that the definition of V here does not reduce to that of V defined earlier where the 
eigenvalues are conventionally ordered oppositely. For example in the MSSM if we had 
required that the CP-even Higgs matrix were diagonalised as VM 2 V^ = diag(m^, m 2 H ) 
then V would have to include a re-ordering of the mass eigenstates, and would have 
been given by 

V = [ ~ Sa y ) (3.13) 

where a is the angle defined earlier in the MSSM. It is to this form that our generalised 
V must reduce in the MSSM limit. 

As in the MSSM we observe that if we do not diagonalise M 2 but instead rotate to 
a basis where the second Higgs doublet does not have any VEV, then the 11 element 
of the matrix gives a useful upper bound on the lightest CP-even Higgs mass 0. This 
rotation is given by WM 2 W^ = M' 2 where 



W 



I Cp Sp 

sp c p j (3.14) 
V 1 



and the bound is given by 



m 2 hi < M[\ = M\ cos 2 2p + AV sin 2 2/3 (3.15) 



The reason why this bound is useful is simply that it has no m' A 2 dependence. For 
example a less useful bound is, 

m 2 hi < Mg = (Mf - AV) sin 2 2(3 + m' A 2 (3.16) 

Starting from the primed basis the CP-even Higgs matrix is diagonalised by 

XM' 2 Xt = diag(r<, ml 2l ml z ) (3.17) 

where X = VWl 

We shall define the relative couplings Ri = Rzzhi as the ZZhi coupling in units 
of the standard model ZZh coupling, and similarly we shall define a ZhiAj coupling 
factor RzhiAj- For example Rzzk^ is a generalisation of sin(/3 — a) and the Rzh^Ai are 
generalisations of cos(/3 — a) in the MSSM. In our notation we find, using the results 
of Ellis et al @, 

Ri = cos (3Vn + sin /3V i2 

= (VW% = X a . (3.18) 
The ZhiAj coupling factorises into a CP-even factor Si and a CP-odd factor Pj|| 

iW; = S^i ( 3 - 19 ) 

where 

= —Vn sin /3 + Vi 2 cos (3 

= (WV%=X i2 (3.20) 

and 

P 3 = (V A U) u+1)2 cos(3 + (V A U) u+1)1 sm(3 (3.21) 

Eq. ( p.21|) implies the simple intuitive results 

Pi = cos 7, P 2 = -sin7 (3.22) 

where 7 defined in Eq. (|3.8| ) is the angle which controls the amount of singlet mixing 
in the CP-odd sector. 
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4 Implementing the Bound on h\ 



In this section we shall consider the absolute upper bound on the mass of the lightest 
CP-even state in this model, using Eq.( [3.15 ) plus radiative corrections. For notational 



ease we shall first define 



Thus the bound is simply 



A = M[j. (4.1) 



< < A 2 (4.2) 

Clearly A 2 is a function of tan/? and A, A 2 (tan/3, A), and to find the absolute upper 
bound we must maximise this function so that 

m 2 hl < A 2 max (4.3) 

where A 2 max is the maximum value of A 2 (tan/3 and A), 

AL, = max (A 2 (tan /3, A)). (4.4) 

Our task in this section is therefore to find A 2 max in the presence of radiative correc- 
tions. 

It is well known that radiative corrections, which we have so far ignored, drastically 
affect the bound 0. In order to deal with these radiative corrections many techniques 
have been proposed 0. Here we shall follow the method proposed in ref. ||, which 
we shall briefly review. According to this method a scale M susy is defined by M^ usy = 
(m~ + m 2 2 )/2 where m^. are the stop mass eigenvalues, and the couplings h t and hb 
found at this scale. Here the squarks are integrated out, leaving an effective Higgs 
potential involving Yukawa couplings with boundary conditions which differ from the 
tree- level ones by some finite corrections ||. Then the Yukawa couplings (including 
h t ) are run from M susy down to m t . Below m t the Yukawa couplings are approximately 
constant, and so the Higgs potential may be minimised at this scale. The effects of 
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the RG running have been estimated analytically, leading to fully analytic results for 
the Higgs masses which agree well with more elaborate methods ||. 

For A = (as in the MSSM) the bound is obviously largest for cos 2 (2/3) = 1. 
Given the restricted range 7r/4 < (3 < 7r/2 this implies that the bound is maximised 
for (3 = 7r/2, cos 2/3 = —1, corresponding to tan/3 = oo (or in practice its maximum 
allowed value). However, for sufficiently large A the tree-level bound will be maximised 
for tan (3 equal to its minimum value. In order to obtain an absolute upper bound, we 
must thus derive an upper limit on A, which will turn out to be a function of tan/3. 
This can be done by demanding that all of the Yukawa couplings remain perturbative 
up to the GUT scale |, |, |. 

We have used the newly discovered mass of the top quark || and the two loop 
SUSY RG equations in the NMSSM |10| to derive an upper bound on A (defined at m t , 



and calculated by requiring all Yukawa couplings ht, lit, and A to remain perturbative 
up to a scale of 10 16 GeV), which we show as a function of tan f3 in Figure 1. Results are 
shown for various values of m t and a^[Mz\ and it is likely that the errors from these 
two measurements will greatly dominate any other uncertainties of our calculation. 
The main features of this graph are that A has an upper bound of around 0.5 to 0.7 
for intermediate tan/3, while for large (small) tan/3 the triviality constraints on hb 
(h t ) force this upper bound very rapidly to zero. 

In Figure 1, we have assumed that the effect of any other Yukawa couplings in 
the singlet sector of the superpotential is negligible, but in fact such couplings can 
have a significant impact, as shown in Figure 2, where we fix tan/3 and explore 
the dependence of the upper bound on A as a function of the coupling k defined in 
Eq. (|3.3| ) for a model with only one singlet. The most important feature of this figure 
is that it clearly displays how adding extra Yukawa couplings will always reduce the 
bound on A, since these will always appear with the same sign in the RG equations. 
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The reason why A falls so rapidly to zero above a certain value of k is that here k is 
approaching its triviality limit, while for small k its impact on the upper bound of A 
is negligible. 

In Figures 3a to 3c, we plot the bound A as a function of tan f3, M susy , and m t 
respectively For each of these figures, we set the squark soft masses squared run, 
m^, and m 2 B to be equal, and we allow the trilinear squark-Higgs coupling A t , A and, 
except in Figure 3a, tan (3 to take on the values which maximise the bound. In Figures 
3a, 3c we select M susy as defined previously to be ITeV, and in Figures 3a and 3b we 
select m t =175GeV. We display bounds for both the MSSM (lower dashed line) and its 
extension with a singlet. Note that for extreme values of tan/3 the maximum allowed 
value of A falls rapidly to zero, and so the two bounds are then the same, while the 
same happens for very large m t since here h t is again close to triviality. Thus the 
large value of m t preferred by CDF [[| means that the two bounds are now within 
around 8GeV. However, it is important to realise that one of the main problems for 
LEP2 will be that for the MSSM small tan f3 always gives a light CP-even Higgs, while 
this need no longer be true in a non-minimal extension, thus seriously damaging the 
prospects for particle searches. 

5 LEP2 Search Strategies 

5.1 Z — > Zhi and General Bounds on CP-even Higgses 

Recently upper bounds on all three neutral CP-even Higgs scalar s in the general 
extension of the MSSM with a single gauge singlet superfield have been derived [0. 
The basic observation which follows from Eq. ( |3.17| ) and the definition M{ 2 = A 2 is: [] 

A 2 = I^m^+lX+XaX 

=_Rlm 2 hl + RW h2 + Rlm 2 h3 (5.1) 

1 Note that X is real if the Higgs sector conserves CP. 
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where we have used the fact that Ri = Xn. Eq.fl5.1|) together with Eq.( |3.12| ) clearly 
implies 

A 2 >(Rl + Rl + Rl)ml (5.2) 

and given that 

R\ + R\ + R\ = 1 (5.3) 
by the unitarity of X, we find the usual bound in Eq. ( |3.15| ), 



m 2 hl < A 2 (5.4) 
Using a similar argument, a bound on the h 2 mass may be obtained from Eq.(5~T) : 



9 A 2 - -R 2 m? . . 

< £ ^ir (5 - 5) 

m| is given from Eq. (|5.1| ) by 



2 iv *H"°hi 1L 2'"'h2 /r c \ 



from which we can extract a bound 



A 2 - (i? 2 + i? 2 )m 2 i 
1 — R\ — R\ 

As we shall see, Eqs. (|5.5| ) and (|5.7| ) are useful upper bounds on and m h3 in terms 
of m\ , i?2 and A 2 = M( 2 (which is just the bound on ml). These bounds tightly 
constrain the spectrum in terms of the couplings and hence will allow us to study the 
reach of colliders when the mixing parameters take on arbitrary values. 

The above general bounds are useful since hi may be light but very weakly coupled. 
In the MSSM this can happen when sin(/3 — a) becomes small. In the present model 
Ri can also be small because hi may contain a large singlet component. In this case 
it does not matter how light hi is since its couplings may be so weak that it can 
never be detected. At first sight this seems to be a catastrophe for LEP2 where in 
the MSSM hi is usually the easiest state to find. However, if R\ pa we may simply 
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ignore h\ and concentrate on hi which then becomes the lightest physically coupled 
CP-even state. Moreover it can be seen from Eq.( |5.5| ) that if R\ — > then hi must 
satisfy 

m 2 h2 < A 2 (5.8) 

which is just the bound on m^. In other words, if the lightest CP-even state is 
essentially invisible, then the second lightest CP-even state must be quite light, and 
may be within the reach of LEP2. Similarly if both R\ and R2 are ~ then Eq.( |5.6| ) 
implies that 

m 2 h3 = A 2 (5.9) 

The above restrictions on the Higgs masses and couplings imply that a collider of 
a given energy and integrated luminosity will be able to exclude values of A below 
which it is impossible for all the Higgses to have simultaneously escaped detection. 

We shall discuss this further in Section 5.3. 

5.2 Exclusion Plots in the R 2 — rrih Plane at LEP 

In this sub-section we we shall be concerned with exclusion plots in the R 2 — m^ 
plane such as those that have already been obtained at LEP1 and whose possible 
form at LEP2 has been predicted, where R 2 = sin 2 (/3 — a) in the MSSM. LEP1 has 



not discovered a CP-even Higgs boson [11], and this non-observation has enabled a 



95% CL exclusion limit to be extracted on the value of the square of the relative 
ZZh coupling R 2 , as a function of the CP-even Higgs mass m^ extending out to 



mh = 65 GeV|12|. Similar exclusion plots at LEP2 will extend the Higgs mass range 
to m h > 65 GeV leading to possible exclusion curves whose precise shape will depend 
on the energy and integrated luminosity of LEP2JTBJ. Three different sets of LEP2 



machine parameters have been considered corresponding to energies and integrated 
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luminosities per experiment of [13 



y/s = 175 GeV, J C = 150 pb" 1 
= 192 GeV, J C = 150 pb" 1 

= 205 GeV, J C = 300 pb" 1 . (5.10) 



The existing LEP1 [12| and anticipated LEP2 JT3[ exclusion plots are combined in 
Fig. 4a, for the three different LEP2 scenarios, and assuming a 100% h\ branching 
fraction into bb. The very steep rise of the contours means that, if a standard model- 
like Higgs boson of a certain fixed mass can be produced at LEP (corresponding to 
R 2 = 1) then reducing R 2 ZZh has little effect on its visibility until R ZZh becomes 
quite small. The relatively flat LEP1 regions of the contours correspond to minimum 
threshold values of R ZZ h below which nothing can be seen for any value of the mass. 

The above exclusion plots in Fig. 4a are the result of a sophisticated Monte Carlo 
simulation of the LEP detectors. It is interesting to observe that these contours 
roughly correspond to the value of the coupling R 2 which would yield 50 Higgs events 
of a given mass. The plot of contours of 50 h events as a function of w,h and R 2 
for the three different LEP2 machine scenarios is shown for comparison in Fig.4b. 
Note that this simple parameterisation of the exclusion plots of Fig. 4a fails badly for 
rrih « M z , and also the actual exclusion limits from LEP1 are significantly stronger 
than our LEP2 approximation for small Higgs masses. Nevertheless Fig. 4b provides 
a useful caricature of the exclusion plots in Fig.4a. 

Finally note that the exclusion plots in Fig. 4 may be interpreted for each of the 
three CP-even Higgs bosons in this model taken separately. In other words the 
disallowed region to the upper and left of the lines will exclude values of Ri 2 for a 
given m/u, and similarly for the other two CP-even Higgs bosons. 
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5.3 Excluded Values of A 



Although there are a large number of parameters in this model, and three CP-even 
Higgs bosons, it proves possible to exclude values of A smaller than a certain amount, 
where A is defined as above to be the upper bound on the lightest CP-even Higgs 
boson, and is to be regarded as a function of A, tan/3 and all the other parameters 
which enter the calculation of radiative corrections such as the top mass, the squark 
masses, and so on. If LEP2 does not discover any CP-even Higgs bosons then exclusion 
plots such as those in Fig. 4a may be produced. We now show that such exclusion 
plots may be used to place an excluded lower limit on the value of A in this model. 
If the excluded lower limit on A reaches the theoretical upper limit of about 146 
GeV (dependent on the top mass and SUSY spectrum as discussed above) then the 
model is excluded. For example, a linear collider of energy 300 GeV should be able 
to exclude the model @. We now describe how LEP may be used to exclude values 
of A in this model. 

Clearly a specified value of A is consistent with many sets of values of the pa- 
rameters rrihi-, Ri subject to the bounds discussed in section 5.1. In general some of 
these sets of m^., Ri will lead to one or more CP-even Higgs boson in the excluded 
(upper left) part of Fig.4a and some sets will lead to Higgs bosons only in the allowed 
region. This means that, for a given A, one or more CP-even Higgs bosons may or 
may not be discovered at LEP2, depending on the values of the other (complicated 
unknown) parameters in the model. According to our discussion in Section 5.1, it is 
clear that as A is reduced, more and more of the allowed sets of m^, Ri will move 
into the excluded part of Fig.4a. As A is reduced below some critical value all the 
sets of values of rrih^Ri will eventually fall into the excluded region of Fig. 4a. This 
critical value of A is the maximum value of A excluded by LEP. This implies that 
the maximum excluded value of A is determined by whether the "worst case" {i.e. 
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hardest to see experimentally) values of m^, Ri consistent with this value of A lie 
within the allowed region or not. 

In order to determine the "worst case" parameters we first fix A at some specified 
value (less than its maximum as determined in section 4) and then we scan over values 
of from zero up to A. For each value of we set R\ equal to the maximum 
allowed value as shown in Fig.4a (or 1 if this would require Ri > 1). For A, 
and R\ fixed as above, we then scan over rrih 2 from up to its upper bound which 
is now fixed according to Eq. (|5.5| ). For each value of rrih 2 we set R<i equal to its 
maximum allowed value as before. For A, m^, m^, Ri and R2 fixed as above, 
and -R3 are now completely specified by Eqs.( |5.3|) and (|5~6l). If R 3 is larger than its 
maximum allowed value (according to Fig.4b), then the "worst case" with the given 
values of m hl and m h2 is excluded. If all the "worst case" points in the m hl , m h2 scan 
turn out to be excluded then we conclude that this value of A is excluded by LEP. 

According to Eq. (|3.15| ), the value of A 2 = M[ 2 is a function of A and tan/3, 
plus parameters which enter in the radiative corrections. For a given value of A, and 
given radiative correction parameters, the exclusion limit on A may be interpreted 
as an exclusion limit on tan/3 in this model, independently of m' A . Thus the present 
analysis will yield a horizontal exclusion line in the m' A — tan f3 plane, as we shall see 
in Section 6. 

If the observability criterion above is replaced by the approximate exclusion limit 
obtained from 50 higgs events regardless of branching fractions, as in Fig.4b, then the 
above algorithm can be solved analytically. In this simple case, the "worst case" for 
LEP2 to discover would be one where the largest of the three cross-sections Oi was 
minimised, where 

Oi = o(e + e~ -> Zhi) = cr SM (e + e~ -> ZK)\ mh=mh% R 2 (5.11) 

It is not hard to show that this will occur when 0\ = 02 = 03, and hence when 
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J2i ( ?SM( m h l )R'i is minimised. Using the constraints in Eqs.flO]) and (|5.3|), together 
with the analytical form for the mass dependence of the tree-level cross-section, 

°SM{m hi ) ~ X 1/2 (X + 12sm|) (5.12) 

where A = (s + (Mz + rrih) 2 )(s — (Mz — ink) 2 ), it is then straightforward to prove that 
this will always be minimised when all three masses are degenerate and R 2 = 1/3. 

This simple result leads to the contour plot of A = \J in the integrated 
luminosity-energy plane shown in Fig. 5. Along each of the contours the energy 
and integrated luminosity correspond to at least one CP-even Higgs boson being pro- 
duced with a yield of 50 events, for any allowed choice of and mh 2 . Thus this 
contour plot is an estimate of the values of A which may be excluded for different 
machine parameters. For example, LEP2 with an energy of 205 GeV will place a 
limit A > 100 GeV, depending on the integrated luminosity. For this energy there is 
a rapid increase in the reach of A from 10-100 GeV as the luminosity is increased from 
100-400 pb -1 , followed by a very slow increase if the luminosity is increased beyond 
this. 

To summarise, the worst case parameters are when all the three CP-even Higgs 

bosons have a mass equal to = A and equally suppressed couplings R 2 = 1/3. In 

fact since in this case there are three Higgs bosons to discover the cross-section will 

be three times larger than for each Higgs boson taken separately. Nevertheless, it is 

clear that for sets of parameters when the Higgs boson masses are not degenerate and 

the Higgs bosons must be considered separately the worst case will never be worse 

than that just described. The excluded value of A is therefore simply determined by 

the following rule of thumb: consider a single Higgs boson with a coupling equal to 

R 2 = 1/3, and find its maximum excluded mass for a given set of machine parameters, 

then equate this mass with the maximun excluded value of A. 

2 Note that this is total luminosity, not luminosity per experiment as elsewhere in this paper. 
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According to this rule of thumb, LEP1 already places a limit on A of A > 59 GeV, 
which is just equal to the mass limit for a CP-even Higgs boson with its ZZh coupling 
suppressed by R 2 = 1/3. Note that because of the steep rises of the function plotted 
in Fig. 4a, the excluded values of A are not far from the present excluded value of 
SM Higgs boson mass of about 65 GeV. Similarly we find that LEP2 will yield the 
exclusion limits A > A min = 81,93,105 GeV, for the three sets of LEP2 machine 
parameters in Eq. (|5.10|) , respectively. 

The above rule of thumb is easily extended to the more general case of n singlets 
being added to the MSSM, i.e. the (M+n)SSM. In such a model the worst case will 
correspond to 2 + n CP-even Higgs bosons each having a mass equal to A, and each 
having a coupling suppression relative to the SM of R 2 = 1/(2 + n). The excluded 
value of A is determined from the exclusion plot of a Higgs boson with a coupling 
equal to R 2 = 1/(2 + n). For example for n = 1, 2, 3 we find that R 2 = 1/3, 1/4, 1/5 
corresponding to the LEP1 excluded values given from Fig. 4a of about A > 59, 58, 57.5 
GeV, respectively. We find it remarkable that such strong limits on A can be placed 
on models with several singlets. 

We emphasise that this argument is only approximate, and may become unreliable 
when more realistic Higgs exclusion data as in Fig.4a are used, rather than the simple 
analytic calculations of how large R 2 may become as a function of m^. (as in Fig4b). 

Finally note that if Higgs bosons are discovered then this does not mean, in 
the context of this model, that A is small. One may have a large A (as large as 
146 GeV ) and still be lucky enough to find that there is a visible CP-even Higgs 
boson corresponding to some fortunate values of parameters. All we have shown is 
how the non-observation of Higgs bosons enables a firm lower limit to be placed on A, 
corresponding to a "worst case" situation. In general, for a given A, the true situation 
will be easier than this, leading to the possibility of Higgs discovery at LEP2 even for 
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very large A. 



5.4 Z -f hA 

In section 3 it was seen that the ZhiAj couplings factorise into a factor from the 
CP-odd matrix multiplied by a factor from the CP-even matrix || given by 

RzhiA-, = S{ cos 7 

RzhiA 2 = -5i£rin7. (5.13) 



Using equation (|3.17| ) we find 



Mg = X 2 2 m 2 hi +X 2 22 m 2 h2 +X 2 2 m 2 h3 

= Sfml + S 2 2 m 2 h2 + S 2 m 2 h3 (5.14) 

where we have used the fact that Si = X i2 . Again 

Sl + S 2 2 +Sl = l (5.15) 

from the unitarity of X. M 22 just corresponds to the less useful bound in Eq.( |3.16|) . 
Eq.( |5.14|) is of course virtually identical to Eq. (|5.1|) , and hence we may immmediately 



write down analagous bounds which relate the masses of the lightest CP-even states 
to their ZhA couplings : 



m 2 hl < M' 2 \ 

2 M 22 — S\m\ 

mi < 



h 2 — i c2 



1-S 2 

, M' 2 \ - (Sf + S 2 )m 2 hi 

m h3 < _ 2 _ 2 (5.16) 

It is trivial to prove that 

m!\ = m 2 Al cos 2 7 + m 2 A2 sin 2 7 (5-17) 
which is similar to Eqs.( |5.1| ) and ( |5.14j ). Eq.( |5.17| ) implies the bound 



m A 2 < m' A 2 (5.18) 
20 



so that m' A is just the upper bound on the lightest CP-odd Higgs boson in this 
model. When A± is weakly coupled (cos7 ~ 0), Eq.( [5.17| ) implies that mA 2 rs rn' A 2 . 
Thus if the lighter CP-odd state is weakly coupled, corresponding to its having a large 
singlet component, then the heavier CP-odd state (which essentially corresponds to 
the MSSM CP-odd state) sits at the bound, in complete analogy with Eq.( |5.6|) . 

We now discuss how LEP limits on the non-observation of Z — > hA may be used 
to exclude values of m' A and M 22 smaller than a certain amount. The following 
discussion obviously parallels that of Section 5.3. As before, the "worst case" for 
Z — ► hA production will provide the best exclusion limit on m' A 2 and M 22 . In order 
to determine this worst case we simply repeat the algorithm in Section 5.3, but with 
the Ri couplings replaced by the Si couplings in Eq.( |5.16| ). However in the present 
case we need to generalise this procedure to include the CP-odd parameters m Al and 
7. Using the exclusion criterion of requiring 50 hA events for exclusion, it is readily 
seen that such a procedure enables the maximum excluded values of m! A and M 22 to 
be obtained. 

Similarly to the case of Section 5.3, we find the worst case to correspond to the 
CP-odd parameters m Al = m A2 = m' A and 7 = 7r/4, and the CP-even parameters 
Si = S2 = S3 = l/\/3 and = = = M 22 . This corresponds to RzhiA 2 = 
1/6 in each case. The excluded values of m' A and M 22 may be found by using the 
approximation of 50 hA events with h having a mass equal to M 22 , A having a mass 
equal to m' A , and the coupling Rzka 2 = 1/6, similarly to the Z — > Zh case discussed 
above. 

As before, the above rule of thumb is easily extended to the more general case 
of n singlets being added to the MSSM, i.e. the (M+n)SSM. In such a model the 
worst case of hA production will correspond to n + 2 CP-even Higgs bosons each 
having a mass equal to J M 22l and each having a coupling factor of S 2 = 1/(2 + n), 
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plus n+l CP-odd Higgs bosons each having a mass equal to m' A , and each having 
a coupling factor of P 2 = 1/(1 + n). As before, when the states are non-degenerate, 
there must be at least one pair hi and Aj which have a cross-section larger than this. 
For example for n — 1, 2, 3, the couplings are RzhA = 1/6, 1/12, 1/20. 

We expect this simple approximation to be rather more robust for the Z — > hA 
process than for Z — > Zh, because of the more rapid fall-off of the cross-section with 
increasing masses for the hA process than for Zh. Furthermore, the larger number of 
parameters here make a more elaborate analysis, such as we shall later perform for 
Z — > Zh, more difficult. Thus we shall always use this approximate method when we 
present our results for Z — > hA. 

We note that of course m' A and M l2 2 are not independent parameters. Using 



Eq. ( |3.16D the value of may be related to m' A and, for a given A and tan (3, M, 



a 
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may be eliminated. Thus, for a fixed value of A, we shall present exclusion plots in 
the m' A — tan (3 plane from the non-observation of Z — > hA. 

5.5 Charged Higgs Detection 

As discussed above, the Z — > hA and Z — > Zh searches are less powerful in extended 
models because the cross-section can be greatly reduced. We now turn to the last 
feasible search at LEP2, that for charged Higgs production. One feature of the model 
which was discussed earlier was that the charged Higgs can be lighter than in the 
MSSM when singlets are included; furthermore, its couplings cannot be suppressed 
by singlet mixing. Hence the charged Higgs signal, which in the MSSM is completely 
dominated by Z — > hA, is now far more important. 

Charged Higgs discovery is complicated, since the rate is strongly dependent not 
only on the Higgs mass, which must not be too close to the W or Z mass to allow 
elimination of background, but also on its branching fractions to cs and rv which are 
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determined by tan j3 |14| . In addition, any LEP2 discovery region may be dominated 
by study of top quark decays at the TeVatron fl5| . For the purposes of this paper 
we shall adopt the rather optimistic view that LEP2 will have sufficient luminosity 
so that the kinematic limit may be approached. 

6 Exclusion Limits in the Plane 

In this section we shall present our results as excluded regions in the m' A — tan f3 
plane, which is familiar from similar studies in the MSSM, and should simplify the 
comparison of the (M+1)SSM with the MSSM. Let us first briefly summarise the 
LEP2 search strategies we have so far introduced in Section 5. 

For the process Z — > hiAj (discussed in section 5.4) we have shown how LEP2 
can be used to place exclusion limits on m' A and M^. For a fixed value of A (and 
radiative correction parameters) these limits may be interpreted as excluded regions 
in the m' A — tan (3 plane. We show both the worst case mixing scenario and the 
simple scenario where there is no singlet mixing. These are calculated using the 
simple approximation that 50 Z — > hA events will be sufficient to ensure a discovery. 

Similarly, for the processes Z — > Zhi (discussed in sections 5.1-5.3) we have shown 
how LEP2 can be used to place an exclusion limit on the value of A. For a fixed value 
of A (and radiative correction parameters) this can be interpreted as an exclusion 
limit on tan/3, independently of m' A (i.e. a horizontal exclusion line in the m^-tan/3 
plane.) However, as we shall see, the resulting excluded region is not very large, so 
we shall resort to a more powerful technique as discussed below. 

This new technique exploits the fact that the upper 2x2 block of the CP-even 
mass squared matrix in Eq.( |3.10| ) is completely specified (for fixed A) in the m^-tan/3 
plane. However, unlike the MSSM, the CP-even spectrum is not completely specified 
since it depends on three remaining unknown real parameters associated with singlet 
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mixing (i.e. the dots in Eq.( |3.10D .) Each choice of these unconstrained terms then 



completely specifies the parameters rrih^Ri, and we can test to see if the resulting 
Higgs spectrum is excluded or allowed. We then scan over all possible choices (which 
we parametrise as , R\ and one other mixing angle) and if the resulting spectrum 
can always be excluded by LEP2, then we conclude that this point in the m' A -tw.j3 
plane (for fixed A) can be excluded. 

Since we have presented two separate algorithms for generating Z — > Zhi contours, 
it is worthwhile here mentioning some of the advantages and disadvantages of the 
different techniques. The method of scanning over the whole of parameter space is 
naturally better than the algorithm based on the exclusion limit of A, in the sense 
that it gives a better reach in the plane. This is clear from the fact that we are 
considering a general 3x3 mass matrix which has the upper 2x2 block given by our 
position in the m' A — tan/5 plane, while the bound algorithm in section 5.3 allows 
any mass matrix which has the correct 11 component in a particular basis (which 
is why the exclusion contours derived in this way do not have any m A dependence). 
On the other hand, the scanning technique is very CPU-intensive, and would rapidly 
become even more so if we allowed extra singlet states. In addition, it is hard to 
understand these results analytically, whereas the simple arguments in section 5.3 
based around Eq. (|5.1|) are much more straightforward. Hence both techniques are 
worth considering. 

6.1 A = 

In figures 6-8 we consider the impact of mixing only, with A set equal to zero. This 
may not be as unreasonably optimistic as it sounds, since recent GUT scale analyses 



|lq , [10]] have concluded that very small A is preferred with universal soft parameters. 
In each case we show the charged Higgs kinematic limit as a dot-dashed line. The 
two Z — > hA exclusion contours are represented by dashed lines and drawn with the 
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assumption of worst case mixing (left) and of no singlet mixing (right). The three 
Z — > Zh exclusion contours are represented by solid lines, and shown for the case of 
no singlet mixing (uppermost line), the full scan over parameter space (middle line), 
and for the bound algorithm (lower horizontal line) discussed in detail in section 5.3. 
The no singlet mixing lines with A = simply correspond to the MSSM. 

Figure 6 shows the exclusion contours for <Js = 175GeV, and integrated luminosity 
of 150pb _1 per experiment. Figure 7 is identical but with ^/s = 192GeV, while Figure 
8 has y/s = 205GeV and integrated luminosity of 300pb _1 per experiment. All these 
figures have m t = 175GeV, degenerate squarks at ITeV, and a 3 (M z ) = 0.12. 

In each of these figures it is clear that the charged Higgs kinematic limit discovery 
line (which, as mentioned above, is rather over-optimistic as a discovery process, but 
is nevertheless indicative of where this process will become important) is completely 
superseded by the Z — > hA line, although to a rather less extreme degree than in the 
MSSM. The inclusion of singlets reduces the excluded region for both Z —>■ hA and 
Z — > Zh quite substantially relative to the MSSM, but still leaves reasonably large 
areas covered. Simply using the bound algorithm however gives rather poorer reach. 

6.2 A > 

To show how large A affects our results, in Figure 9 we show a figure with yfs = 
192GeV, and integrated luminosity of 150pb _1 per experiment, but with A = 0.5. 
Here the charged line is becoming quite competitive with the Z — > hA line, and so 
a charged Higgs search may well be the most practical one at LEP2. The Z — > hA 
contours are hardly changed from those in the A = case except for very small tan/?. 

The most significant impact is however on the Z — > Zh lines. Here even the no 
mixing scenario has very little reach in the rriA — tan/5 plane because, as is clear from 
Eq.( p.l5|) , small tan/3 no longer implies a light CP-even state. With singlet mixing, 
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no part of the plane can be covered at all. 



7 Conclusion 

The addition of extra singlets to the MSSM greatly complicates the model and renders 
far harder the main experimental searches. However it is possible to exclude this 
model by placing an experimental lower bound A min on the value of the parameter A 
which is just the upper bound on the lightest CP-even Higgs boson mass. In the limit 
that the lightest CP-even Higgs boson is very weakly coupled, this bound applies to 
the second lightest CP-even Higggs boson. Thus for a fixed value of A the entire 
CP-even Higgs boson spectrum is constrained, resulting in experimentally excluded 
values of A associated with the "worst case" scenarios discussed in section 5.3. 

We have seen that LEP1 already finds A > A min = 59 GeV, and LEP2 will 
set limits of A min ks 81,93,105 GeV for three different levels of operation. The 
theoretical upper bound on A is A max 146 GeV, depending on the details of the 
squark spectrum and on the top mass. Once A min becomes greater than or equal to 
A max then the model will be excluded. We have generalised this procedure to the 
case of an arbitrary number of extra singlets. Thus our first conclusion is that it will 
be possible to exclude a version of the MSSM containing additional singlets. 

The effects of the additional singlet can be thought of as reducing the Higgs sector 
of the model to the MSSM with two additional complicating factors : extra singlet 
states which can mix in an arbitrary way with the usual neutral Higgs states altering 
the masses and diluting the couplings of the mass eigenstates; and an extra Higgs 
sector quartic coupling A which changes the mass matrices even in the absence of 
singlet mixing. 

We have systematically studied the case of singlet mixing with A = 0, with the 
primary conclusion that, while the inclusion of a singlet can significantly complicate 
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matters for searches at LEP2, it is still possible to cover a significant amount of the 
m' A — tan j3 plane using the usual Z — > Zh and Z — * hA searches. The main difference 
in strategy between the searches in the MSSM and in models with singlets is that 
more luminosity is needed to cover the same area of the plane, since states can be 
more weakly coupled than in the MSSM. A second difference is that, because singlets 
cannot mix with charged states, the charged Higgs signal cannot be degraded and so 
is more important than in the MSSM where it is completely dominated by Z — > hA. 

The effect on non-zero A is more troublesome. For large A the mass bounds in 
the MSSM are markedly increased, particularly at small tan (3 which is of course the 
region where LEP2 normally has the best reach. This largely wrecks the prospect of 
excluding large parts of the m' A — tan (5 plane through the process Z — > Zh, and has 
some impact on Z — > hA. However, we note that large A also reduces the charged 
Higgs mass substantially, making its discovery easier. 

We conclude on a positive note by pointing out that for much of this paper we 
have assumed worst case mixing scenarios, and considered how it might be possible 
to rule out regions of parameter space in a consistent way for arbitrary parameters 
from the singlet sector. For the MSSM, specifying a point in the tua — tan f3 plane 
completely specifies all the masses and mixings, and so points outside the discovery 
contours cannot be discovered; however with singlets it is possible that, if we are 
lucky enough, any point in the m' A — tan (5 plane could lead to a discovery, since it is 
possible to construct singlet mixing parameters which will reduce the mass of a heavy 
CP-even state while leaving its coupling reasonably large. Hence we can argue that 
despite the rather negative impact of the mixing on the LEP2 exclusion contours, 
a SUSY Higgs discovery may be no less likely in a model with singlets than in the 
MSSM. 
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8 Figure Captions 



Figure 1: Maximum value of A(m t ) as a function of tan (3, shown for various values 
of m t and a 3 (M^) and derived from requiring perturbativity of all Yukawa couplings 
up to a scale 10 16 GeV. 

Figure 2: Maximum value of \{m t ) as a function of k{m t ), shown for various val- 
ues of tan/5, with m t = 175GeV and a 3 (Mz) = 0.12 and derived from requiring 
perturbativity of all Yukawa couplings up to a scale 10 16 GeV. 

Figure 3a: A max , the bound on the lightest Higgs in the MSSM with (upper solid 
line) and without (lower dashed line) added singlets, shown as a function of tan/3 
with M susy = ITeV, m t = 175GeV, and A and A t set so as to maximise the bound. 

Figure 3b: A rnax , the bound on the lightest Higgs in the MSSM with (upper solid 
line) and without (lower dashed line) added singlets, shown as a function of M susy 
with m t = 175GeV, and A, A t and tan/5 set so as to maximise the bound. 

Figure 3c: A rnax , the bound on the lightest Higgs in the MSSM with (upper solid 
line) and without (lower dashed line) added singlets, shown as a function of m t with 
M susy = 175GeV, and A, A t and tan/5 set so as to maximise the bound. 

Figure 4a: R 2 versus m h for three different sets of machine parameters at LEP2, 
using a full analysis of the experimental data. For m h ~50GeV the LEP1 bounds are 
tighter and so we use those instead. These are 95% confidence level exclusion contours. 
Points above the line correspond to states which could be seen, those below to those 
which could not be seen. The three lines are for y/s = 175GeV, / C = 150pb _1 
per experiment; \/s = 192GeV, J C — 150pb _1 per experiment; y/s = 205GeV, 
/ C = 300pb _1 per experiment; reading from left to right. 

Figure 4b: R 2 versus rrih for three different sets of machine parameters at LEP2, 
using the simple approximation that 50 CP-even events are sufficient for detection of 
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a state. Points above the line correspond to states which could be seen, those below 
to those which could not be seen. The three lines yfs = 175GeV, / C = 150pb _1 
per experiment; \fs = 192GeV, f C — 150pb _1 per experiment; y^s = 205GeV, 
J C = 300pb _1 per experiment; reading from left to right. 

Figure 5: Contours of values of A min , the value of A which can be excluded, in the 
y/s— luminosity plane, labelled in GeV. This is total luminosity, not luminosity per 
experiment. 

Figure 6: m' A — tan (3 plane with v /i=175GeV, f C — 150pb _1 per experiment, 
A = 0, m t = 175GeV, and squarks degenerate at ITeV. The dot-dashed line is the 
kinematic limit for charged Higgs discovery. The dashed lines represent Z — > hA 
contours with no singlet mixing (right, corresponding to the MSSM) and worst case 
mixing (left). The three solid lines show Z — > Zh contours, derived by assuming no 
singlet mixing (upper, corresponding to the MSSM), the full scan approach (middle) 
and the exclusion limit on A (lower). 

Figure 7: m' A — tan (3 plane with v /s=192GeV, / C = 150pb _1 per experiment, 
A = 0, m t = 175GeV, and squarks degenerate at ITeV. The dot-dashed line is the 
kinematic limit for charged Higgs discovery. The dashed lines represent Z — > hA 
contours with no singlet mixing (right, corresponding to the MSSM) and worst case 
mixing (left). The three solid lines show Z — > Zh contours, derived by assuming no 
singlet mixing (upper, corresponding to the MSSM), the full scan approach (middle) 
and the exclusion limit on A (lower). 

Figure 8: m' A — tan/3 plane with v /s=205GeV, J C — SOOpb" 1 per experiment, 
A = 0, m t = 175GeV, and squarks degenerate at ITeV. The dot-dashed line is the 
kinematic limit for charged Higgs discovery. The dashed lines represent Z — * hA 
contours with no singlet mixing (right, corresponding to the MSSM) and worst case 
mixing (left). The three solid lines show Z — > Zh contours, derived by assuming no 
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singlet mixing (upper, corresponding to the MSSM), the full scan approach (middle) 
and the exclusion limit on A (lower). 

Figure 9: m' A — tan/3 plane with v /s=192GeV, f C — 150pb _1 per experiment, 
A = 0.5, m t = 175GeV, and squarks degenerate at ITeV. The dot-dashed line is the 
kinematic limit for charged Higgs discovery. The dashed lines represent Z — > hA 
contours with no singlet mixing (right) and worst case mixing (left). The solid line 
shows the Z — ► Zh contours derived by assuming no mixing. The full scan approach 
and using the limit on A do not allow any part of the plane to be excluded. 
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